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A potential can have features that do not reflect the dynamics of the system it 
describes but rather arise from the choice of interpolating fields used to define it. 
This is illustrated using a toy model of scattering with two coupled channels. A 
Bethe-Salpeter amplitude is constructed which is a mixture of the waves in the two 
channels. The potential derived from this has a strong repulsive core, which arises 
from the admixture of the closed channel in the wave function and not from the 
dynamics of the model. 

I. INTRODUCTION 

In recent work on two-baryon systems in lattice QCD, members of the HAL-QCD collab- 
oration have extracted baryon-baryon potentials from their simulations [IHZ]- The method 
relies on the construction of a Bethe-Salpeter amplitude or "wave function" - a matrix el- 
ement of two interpolating fields with the quantum numbers of the baryons involved. This 
amplitude is then inserted into a Schrodinger equation in order to deduce the baryon-baryon 
potential that would generate it. More recently, the same idea has been used to extract short- 
distance potentials through the operator product expansion of perturbative QCD [BHTT]. A 
review of the approach can be found in Ref. [12] . 

One interesting feature of the resulting nucleon-nucleon potential is that it possesses a 
repulsive core pQ, very reminiscent of those in many of the potentials traditionally used 
in nuclear physics, such as Argonne vi% [13J. When extended to three flavours of quark 
with SU(3) symmetry, repulsive cores are seen in all channels except the SU(3) singlet [I]. 
The repulsion is particularly strong in the symmetric octet channel, where no long-range 
attraction is seen. 

This pattern is qualitatively similar to that seen in constituent quark models (for reviews, 
see Refs. [HHPo] ). For two octet baryons, the symmetric octet channel is forbidden by the 
Pauli principle at small separations, leading to strong repulsion. Other channels are not 
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forbidden and so, as well as antisymmetrisation, the form of the interaction between the 
quarks is needed, as discussed by Oka [T7]. In quark models based on magnetic gluon 
exchange, short-distance repulsion is found in all channels except the singlet. 

The HAL-QCD approach has been criticised, for example in Refs. [181 HH], because a 
potential is not an observable in either experiments or lattice simulations. In the lattice 
case, the energies of two-baryon states are observables and these can be related, via Luscher's 
formula [20J, to experimental observables, namely, phase shifts. In the HAL-QCD approach, 
this information on the physical phase shifts is encoded in the tails of the Bethe-Salpeter 
amplitudes outside the interaction region. In contrast, the short-distance forms of these 
amplitudes, and hence the deduced potentials are dependent on the choice of interpolationg 
field used to define them. 

In this short note, I use a toy model of scattering with two coupled channels to illustrate 
how a potential constructed in this way can develop features that reflect the choice of 
interpolating fields, and not the actual dynamics. In particular an admixture of a closed 
channel can lead to a strong repulsive core in the derived potential that does not correspond 
to any aspect of the potentials in the original model. 

II. COUPLED-CHANNEL MODEL 

The toy model I use to illustrate issues with interpreting potentials from Bethe-Salpeter 
wave functions is two-body scattering in one dimension with two coupled channels. The 
lower-energy channel can be thought of as corresponding to two nucleons in their ground 
states; the other to one ground-state nucleon and one N* with the same quantum numbers. 
I shall refer to the channels as "NN" and "NN*". The particles in the NN channel interact 
through an attractive square-well potential similar to the long-range attractive force between 
nucleons. For simplicity I take the potential in the NN* channel to be zero. The N* has 
an excitation energy A and so the NN* channel is closed for energies below A. The two 
channels are coupled by a contact interaction of strength g. 

While this model is very much a caricature of real nucleon-nucleon scattering, I believe 
that it retains enough aspects of the real system to illustrate why one should be wary about 
taking seriously short-distance features. 

The wave functions for the relative motion of the particles in the two channels, ipo{x) and 
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ipi(x), satisfy the coupled equations, 

1 d 2 



M dx 2 
'Mdx 2 



+ V(x) - E 
+ A-E 



i) (x)+2g5(x)i) 1 {x) = 0, 
ij)i(x) + 2gS(x)ip (x) = 0, 



(1) 



where 



V{x) 



(2) 



— Vo for \x\ < a 
for |x| > a. 

A general interpolating field for the "nucleon" in this model is just a linear combination 
of the N and N* fields. A Bethe-Salpeter amplitude *f?(x) can be defined for any coupled- 
channel solution by first building a second-quantised state representing the solution and 
then taking the matrix element of a pair of interpolating field operators between this state 
and the vacuum. The result is just a linear combination of the two wave functions, 



^(x) = ip (x) +pipx(x). 



(3) 



Since it is not required for the determination of the potential, I have not specified the overall 
normalisation of this quantity; only the coefficient p of the N* admixture is relevant. 
The (even-parity) solutions of the coupled equations ^ have the forms 



ipi(x) 



Acos(K\x\) + B sm(K\x\) for \x\ < a 

C cos (A; | a; |) + D sm{k\x\) for \x\ > a, 
exp(— a\x\), 



(4) 



where 



k 2 
K 2 
a 2 



ME, 

M(E + V ) 
M(A-E). 



(5) 



Again the overall normalisation is irrelevant, and so I have arbitrarily taken ipi(0) = 1. 

Matching the solutions at the edge of the square well and at the origin, where the 8- 
function acts, leads to the boundary conditions, 



V'o(o-) = i ) o(a+), 
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^o(a-) = ^o( a +)> 
$,(0) = MgfaiO), 

^(0) = MgMO)- (6) 
These give rise to a set of linear equations for the coefficients, A, B, C, D: 

A cos(Ka) + B sm(Ka) = C cos(ka) + D sm(ka), 

—KA sm(Ka) + KB cos(Ka) = —kC sin(ka) + kD cos(fca), 

KB = Mg, 

-a = MgA, (7) 

which can be solved straightforwardly. The resulting wave functions at energy E will be 
denoted 4>o,i(x; E). From these, we can construct the Bethe-Salpeter amplitude *ff(x; E). 

Following the approach outlined in the introduction, I now use *ff(x]E) to define a po- 
tential by demanding that it satisfy the Schrodinger equation, 

1 d 2 



Mdx ^Vss^E)-E 



V{x;E) = Q. (8) 



The resulting potential is 



III. RESULTS 

I present results here for a representative parameter set, MVq = 1 /a 2 , MA = 6/ a 2 and 
M g = Q/a, and I take a = 1 for convenience in plotting the results. Fig. 1 shows the zero- 
energy potential Vbs{x~, 0) for a range of values of the mixing parameter in the interpolating 
field. Note that I have not plotted the potential for x > a, where it is essentially zero, nor 
have I shown the 8 function that is present at the origin. 

A couple of features are worth noting about these potentials. First, the long-range at- 
traction in the NN channel, V(x) = —Vq, can be seen as x approaches 1. This similar to the 
way that the pion-exchange tail can be seen in the potentials extracted by the HAL QCD 
collaboration [TJ. The second feature is the apparent repulsive "hard core" that is present 
in all cases for x ^5 0.5. This comes entirely from the closed NN* channel, which makes a 
large positive contribution to Vbs{x] 0) for x < 1/a because d 2 ipi(x) / dx 2 ~ MA^(i) there. 
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FIG. 1: The zero-energy potential MV BS (x;0) for p = 0.1, 0.2, 0.3, 0.5, 1, 2, 3. 

It has nothing to do with the NN-channel wave function, as demonstrated by its strong de- 
pendence on the choice of interpolating field. In contrast, the long-range attraction (where 
it is not swamped by the core) is independent of the choice of field. 




FIG. 2: The potential MV B s(x; E) for ME = 0, 1, 2, all with p = 1. 

At low energies, the core in this potential is only weakly dependent on the energy E, as 
can be seen in Fig. 2. Like the HAL-QCD potentials, this has the appearance of something 
that could be well approximated by a local potential. However this is just a reflection of 
the fact that the strength and range of the core arise from the closed-channel wave function 
i/ji(x), and so are controlled by g and A (and p of course), so long as E lies well below the 
excitation energy A. 
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IV. CONCLUSION 

The results in the previous section illustrate how a potential can have features that 
arise from the interpolating fields used to define it. In particular they show how coupling 
to closed channels can lead to repulsive cores in potentials derived from Bethe-Salpeter 
amplitudes. Similar repulsive cores seen in the HAL-QCD collaboration's analyses of their 
lattice simulations [IHZ] could therefore be results of their choice of interpolating fields, 
rather than indications of similar dynamical mechanisms to those of the constituent quark 
model. 

Since the same interpolating fields are used in all channels, one might argue that the 
absence of a repulsive core in the SU(3)-singlet channel is evidence against these cores being 
artefacts of the choice of field. However the strength of the core is dependent on the size 
of the coupling to the closed channel and, for baryons in the SU(3) limit this coupling is 
likely to be strongly channel-dependent. In particular, this dependence is to be expected if 
the coupling to excited octet baryon states occurs via intermediate excitation of decuplet 
baryons. Such excitations could be mediated by pion and kaon exchanges at long distances 
or, at shorter ranges, by the same magnetic gluon interactions that are responsible for the 
repulsive cores in the quark model [T7j. This mechanism for generating excited octet baryons 
would not operate in the SU(3)-singlet channel since the relevant intermediate states - octet- 
decuplet or decuplet-decuplet - cannot couple to an overall singlet. It could therefore explain 
the absence of a repulsive core in that channel. 

To summarise: a potential is not an observable. The toy model in this note provides a 
warning that even quite striking features of a potential may not reflect the actual dynamics 
of the system it describes. 
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